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Abstract. An abstract should be given 



In this Letter we further explore the idea, suggested previ- 
ously by Kluzniak and collaborators, that the high frequency 
QPOs may be explained as a resonant coupling between the 
neutron star spin and two epicyclic modes of accretion disk os- 
cillations. We confirm result of Lee, Abramowicz and Kluzniak 
J2004i that the strongest response occurs when the frequency 
difference of the two modes equals either the spin frequency 
(for "slow rotators") or half of it (for "fast rotators"). New 
points discussed in this Letter are: (1) We suggest that the cou- 
pling is gravitational, and due to a non-axially symmetric struc- 
ture of the rotating neutron star. (2) We found that two excited 
modes may be both connected to vertical oscillations of the 
disk, and that strong gravity is not needed to excite the modes. 
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1. INTRODUCTION 

To date, quasi-periodic oscillations (QPOs) have been observed 
in about twenty Low Mass X-ray Binaries (LMXBs) sources 
containing an accreting neutron star. Among these systems, 
the high-frequency QPOs (kHz-QPOs) which mainly show up 
by pairs, denoted by frequencies vi and V2 > vi, possess 
strong similarities in their frequencies, ranging from 300 Hz to 
about 1300 Hz, as well as in their shape (see van der Khs 120001 
for a review). For slow rotators (i.e. with rotation rate v, ^ 300 
Hz), the frequency difference between the two peaks is around 
Av w V, and V2 ~ 3 vi /2 whereas for fast rotators (i.e. v, ^ 600 
Hz), this difference is around Av « v»/2, (van der Klis '2004i. 
For black hole candidates the 3:2 ratio was first noticed by 
Abramowicz and Kluzniak <2001> who also recognized and 
stressed its importance. Now the 3:2 ratio of black hole QPOs 
frequencies is well established (McChntock and RemiUard, 

HOE)- 



Many attempts have been made to explain this phe- 
nomenology. The relativistic precession model introduced by 
Stella & Vietri ( 1998. .199 9 1 makes use of the motion of a sin- 
gle particle in the Kerr-spacetime. However, the peak separa- 
tion is not naturally deduced from their model. More promis- 
ingly, Abramowicz & Kluzniak ( 2 001 J introduced a resonance 
between orbital and epicyclic motion that can account for the 
3:2 ratio around KeiT black holes leading to an estimate of 
their mass and spin. Kluzniak, Abramowicz, Kato et al. (.2004a> 
showed that the twin kHz-QPOs is explained by a non Un- 
ear resonance in the epicyclic motion of the accretion disk. 
Rebusco (2004 1 developed the analytical treatment of these 
oscillations. Bursa, Abramowicz, Karas & Kluzniak J2004> 
suggested a gravitational lens effect exerting a modulation of 
the flux intensity induced by the vertical oscillations of the 
disk while simultaneously oscillating radially. More recently, 
Torok, Abramowicz, Kluzniak & Stuchlik ( 2005i applied this 
resonance to determine the spin of some microquasars. 



It was recognized (Kluzniak et al. I2004al Kluzniak, 
Abramowicz & Lee I2004b> that because in the accreting mil- 
lisecond pulsar SAX J1808.4-3658 the difference in frequen- 
cies of the double peaked QPOs is clearly equal to half of the 
pulsar spin (Wijnands, van der Klis, Homan et al. I2003t . the 
epicyclic resonance must be excited by a coupling of the accre- 
tion disk oscillation modes to the neutron star spin. Numerical 
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simulations by Lee, Abramowicz & Kluzniak J2004t modeUed 
the coupling by an unspecified external forcing of the disk 
(with periodicity equal to that of the spin) and found that reso- 
nant response occurs when the difference between frequencies 
of the two modes equals to one-half of the spin frequency (as 
observed in SAX J 1808.4-3658 and other "fast rotators"), and 
when it equals to the spin frequency (as observed in "slow ro- 
tators" like XTE J1807-294). 

In this Letter we further explore these ideas by showing 
that the desired coupling may be provided by gravity of a suf- 
ficiently non-axially symmetric neutron star. (Discussion of a 
physical plausibility of non-axially symmetric neutron stars 
is beyond the scope of this Letter). In such a case, the ac- 
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cretion disk will experience a gravitational field with dipolar, 
quadrupolar and octupolar moments (m= 1,2,3) that vary peri- 
odically in time. 

2. THE MODEL 

In this section, we describe the main features of the model, 
starting with a simple treatment of the accretion disk, assumed 
to be made of non interacting single particles orbiting in the 
equatorial plane of the star. We thus neglect the hydrodynam- 
ical aspect of the disk like pressure. Particles evolve in a per- 
fectly spherically symmetric gravitational potential until, at the 
origin of time t - 0, a rotating asymmetric part is added to 
the stellar gravitational field. This perturbation is issued from 
an inhomogeneity in the neutron star interior, for instance, as 
explained in the following subsection. 

2.1. Distorted Stellar gravitational field 

We assume that the stellar interior is inhomogeneous and 
anisotropic. In some regions inside the star, clumps of mat- 
ter generate locally a stronger or weaker gravitational potential 
than the average, depending on their density. In order to com- 
pute analytically such kind of gravitational field, we idealize 
this situation by assuming that the star is made of an homo- 
geneous and isotropic matter everywhere (with total mass M, 
and angular speed To this perfect spherically symmetric 
geometry, we add a small mass point, the perturber having a 
mass Mp «: M» located inside the star at a position /?p = 
(rp,<^p = r,Zp) with + < RI- We use cyhndrical co- 
ordinates denoted hy R - (r, (p, z). The origin of the coordinate 
system coincides with the location of the neutron star. A finite 
size inhomogeneity can then be thought as a linear superpo- 
sition of such point masses. The total gravitational potential 
induced by this idealized rotating inhomogeneous star is : 



(D(r, (fi, z,t)^-GMJ\\R\\-GMp/\\R- Rp 



(1) 



where the first term in the right hand side corresponds to 
the unperturbed spherically symmetric gravitational potential 
whereas the second term is induced by the small point like 
inhomogeneity. Using the cylindrical frame of reference, the 
gravitational potential reads : 



^{r,ip,z,t) - -G Mtl yjr^ + z^ -GMpI ^jr^ + z^ x 

2 ' ( ^rl+zH ^lr^+z^) cos(m ^) (2) 



m=0 



where the azimuth in the corotating frame is, ip - ^p - Q.t t and 
the Laplace coefficients b'l^{x) of celestial mechanics are : 
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cos(m 41) 
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dij/ 



(3) 



where (jj), is the Kronecker symbol. The total linear response of 
the disk is then the sum of each perturbation corresponding to 
one particular mode m. Because the perturber is inside the star 
and the disk never reaches the stellar surface, x < 1 and thus the 

1 /2 

Laplace coefficients (x) never diverge. Moreover, because 



of the term cos(m 1/') in the integrand Eq. Q, the value of the 
Laplace coeflicients decreases rapidly with the azimuthal num- 
ber m. As a result, only the low azimuthal modes will influence 
significantly the evolution of the disk. Keeping only the few 
first terms in the expansion is sufficient to achieve reasonable 
accuracy. 

2.2. Equation of motion for a test particle 

All particles evolve in the gravitational field imposed by the 
neutron star. Their equation of motion reads 



r-r(p = gr + Sgr 
2r(p + rip = 6g^ 



(4) 
(5) 
(6) 



The gravitation field of the star M, is denoted by g whereas 
that of the perturber Mp is denoted by 6g and the dot means 
time derivative didt. We are only interested in the vertical mo- 
tion experienced by the test particles in response to the verti- 
cal perturbed gravitational field. We therefore neglect the ra- 
dial and azimuthal perturbations, 6gr - dg^ - 0. We keep 
only the vertical component, Sgz + 0. According to this sim- 
plification, the second equation then integrates immediately. It 
states the conservation of the angular momentum of the parti- 
cle L = m = const. This is the obvious integral of motion 
for this problem. 

Perturbing Eq. (|6j and developing to first order in the per- 
turbation around the equilibrium Keplerian orbit defined by 
(r = ro, 1^ = Qk Z() - 0), the vertical motion reads : 



Z H- 



(7) 



m=0 
+00 

m=0 

where the Laplace coefficient are evaluated at the point (ro, zo), 

= bli'(^f(^p/^[r[^, fik = (GMJriy'^ is the 

Keplerian orbital frequency and Qp = (G Mp/r^y^^. We recog- 
nize a Hill equation (periodic variation of the eigenfrequency 
of the system on the left hand side) with aperiodic driving force 
(on the right hand side). 

2.3. Resonance conditions 

Eq. Q describes an harmonic oscillator with periodically vary- 
ing eigenfrequency which is also excited by a driven force. It 
is well known that some resonances will therefore occurs in 
this system. Namely, we expect three kind of resonances corre- 
sponding to : 

- a corotation resonance at the radius where the angular ve- 
locity of the test particle equals the rotation speed of the 
star. This is only possible for prograde motion. The reso- 
nance condition is Ok = Q,; 

- a driven resonance at the radius where the vertical epicyclic 
frequency equals the frequency of each mode of the grav- 
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itational potential as seen in the locally corotating frame. 
The resonance condition is m |Q, - Qk| = k^; 
- a parametric resonance related to the time-varying vertical 
epicycUc frequency, (Hill equation). The rotation of the star 
induces a sinusoidally variation of the vertical epicyclic fre- 
quency leading to the well known Mathieu's equation for a 
given azimuthal mode m. The resonance condition are : 



m |f2t - Qkl -2 — 
n 



(8) 



where n> 1 is a natural integer Note that the driven resonance 
is a special case of the the parametric resonance for n - 2. 
However, their growth rate differ by the timescale of the am- 
plitude magnification. Driving causes a linear growth in time 
while parametric resonance causes an exponential growth. We 
also rewrite the vertical epicyclic frequency as instead of f2k 
in order to apply the results to a more general case which could 
include magnetic field or general relativistic effects. Indeed, in 
Sec|3]we will apply the aforementioned resonance criteria also 
to the Kerr spacetime for which the degeneracy between orbital 
and vertical epicyclic frequency is lifted (/Cz Qj^). 

3. RESULTS 

3.1. Newtonian disk 

From Eq. (|8}, we can find the radius where each of this reso- 
nance will occur Beginning with the Newtonian potential, it is 
well known that the angular velocity and the vertical epicyclic 
frequencies for a single particle are equal so that Qk = k^,. 
This conclusion remains true for a thin accretion disk hav- 
ing Cs/rQk <K 1 where c^ is the sound speed. Distinguishing 
between the two signs of the absolute value, we get for the 
parametric resonance condition Eq. the following orbital 
rotation rate : 

Qk m 
Qt m + 2/n 



(9) 



As a consequence, the resonances are all located in the fre- 
quency range Qk £ [Q»/3,3f2»]. In Tabled we indicate the 
results for a slow as well as for a fast rotator (respectively 
V, = 300 Hz and v» - 600 Hz). Numerical applications 
are given for a spinning neutron star, showing the first three 
modes m and the first two integers n = 1,2. The pair of high- 
est orbital frequencies for the v, = 300 Hz spinning neutron 
star are vi = 600 Hz and V2 - 900 Hz. The twin peak separa- 
tion frequency is then Av = 300 Hz - v». The vertical motion 
induced by the parametric resonance at that location will ap- 
pear as a modulation in the luminosity of the accretion disk. 
For the v» = 600 Hz spinning neutron star, the highest orbital 
frequencies are 1800 Hz and 1200 Hz. However, due to the 
ISCO, the former one is not observed because it is located in- 
side the ISCO and therefore does not correspond to a stable or- 
bit. Indeed, remind that for a 1 .4 M© neutron star, the maximal 
orbital frequency in the Schwarzschild spacetime at the ISCO is 
Visco - 1571 Hz. Therefore, the first two highest observable fre- 
quencies are vi - 900 Hz and V2 - 1200 Hz. This is confirmed 
in the next subsection where the Newtonian field is replaced 
by the Kerr geometry. The peak separation frequency becomes 



Table 1. Value of the orbital frequencies at the parametric res- 
onance for the first three modes m and with n = 1 , 2, in the case 
of a Newtonian gravitational potential. The results are given for 
a 1.4 Mq neutron star rotating respectively at 300 and 600 Hz. 
The value on the left of the symbol / corresponds to the absolute 
value sign taken to be - and on the right to be -H. 



Mode m 



Orbital frequency v{r, a,) (Hz) 





V, = 600 Hz 


V. = 300 Hz 




n = I n = 2 


n = I n = 2 


1 


-600 / 200 — / 300 


-300/ 100 — / 150 


2 


— /300 1200/400 


— /150 600/200 


3 


1800/360 900/450 


900/ 180 450/225 



then Av = 300 Hz - v,/2. Thus the peak separation for slow 
spinning neutron stars is Av = v* whereas for fast spinning 
neutron star it becomes Av = v,/2. This segregation between 
slow and fast rotating neutron stars is well observed in several 
accreting systems (van der Klis 20041. A similar reason for the 
dichotomy between fast and slow rotators (location of the res- 
onance radius), was previously suggested by Lee, Abramowicz 
& Kluzniak (l2004l. 



3.2. General relativistic disk 

When the inner edge of the accretion disk reaches values of a 
few gravitational radii, the general relativistic effects become 
important. The degeneracy between the frequencies Qk and 
is lifted and will depend on the angular momentum of the 
star The characteristic orbital and vertical epicyclic frequen- 
cies in the accretion disk around a Kerr black hole (or a ro- 
tating neutron star) are the orbital angular velocity, Q(r, a,) - 
GMJ(P'^ + Rl'^a^) where Rg = GMJc^ is the gravita- 
tional radius and the vertical epicyclic frequency, Kj{r,a») - 

Q.{r, fl,) ^1-4 (R^lrfl^ a, -H 3 {R^jrY al The parameter a, 
corresponds to the angular momentum of the star, in ge- 
ometrized units. For a neutron star of mass Mt, moment of 
inertia /» and rotating at the angular velocity Q.,, it is given 
by fl» - c It D.t/G M^. The parametric resonance conditions 
Eq. (|8} splits into two cases, depending on the sign of the ab- 
solute sign : 



Q(r,fl.) + 



2 K^ir, a.) 



(10) 



For a given angular momentum a,, Eq. ilQ\ have to solved for 
the radius r. For a neutron star, we adopt the typical parame- 
ters : 



- mass M« = 1 .4 Mq ; 

- angular velocity v, - 

- moment of inertia /, 



Q./27r = 300/600 Hz; 
= 10^^ kgm^ ; 



The angular momentum is then given by a* - 5.79 x 10"^ Q,. 
Solving Eq. ilO\ for the radius and then deducing the orbital 
frequency at this radius we get the results shown in Table |2 
For the above chosen stellar spin rate we find a, - 0.1/0.2 
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for the slow and fast rotator respectively. Therefore the vertical 
epicyclic frequency remains close to the orbital one k^^ ^ Q^. 
As a consequence, for the vertical resonance, the Newtonian 
approximation mentioned in the previous subsection remains 
valid. Indeed, for a slow rotator, the two highest frequencies 
are vi = 599 Hz and V2 = 898 Hz, therefore Av = 299 Hz 
which is still very close to v» whereas for a fast rotator, v\ - 
899 Hz and V2 = 1198 Hz, therefore Av = 299 Hz which is 
close to v,/2 (the orbit located within the ISCO has been dis- 
carded). Numerical simulations in which the disk was disturbed 

Table 2. Same as Table [l] but in the general relativistic Kerr 
space time. 



Mode m Orbital frequency v(r, a,) (Hz) 





V* = 


600 Hz 


V* = 


300 Hz 




n = 1 


n = 2 


n = 1 


n = l 


1 


/ 200 


— / 300 


— / 100 


— / 150 


2 


— / 300 


1198/400 


— / 150 


599 / 200 


3 


1790 / 360 


899/450 


898/ 180 


450/225 



by an external periodic field confirmed this point of view (Lee, 
Abramowicz & Kluzniak 2004 1. We emphasize the fact that 
these results apply to a rotating asymmetric magnetic field with 
exactly the same resonance conditions Eq. ^ provided that 
the flow is not to far from its Keplerian motion, i.e. a weakly 
magnetized accretion disk with high y6-plasma parameter Note 
also that these results are quiet general and independent of the 
spacetime geometry. General relativity is not required to ac- 
count for the 3:2 ratio. This model therefore also encompasses 
the accreting white dwarfs for which QPOs have been observed 
in the same 3:2 ratio and show strong similarities with X-ray 
binaries, (Warner & Woudt l2005l . 

4. CONCLUSION 

In this Letter, the consequences of a weak rotating asymmet- 
ric gravitational potential perturbation on the evolution of a 
thin accretion disk initially in a stationary axisymmetric state 
have been explored. For gravitational perturbation with multi- 
polar components, the response of the disk is the sum of indi- 
vidual modes as long as it remains in the linear regime. The 
physical processes at hand does not require any general rela- 
tivistic effect. Indeed, the resonances behave identically in the 
Newtonian as well as in the Kerr field. As a consequence, the 
QPO phenomenology is unified in a same picture, whatever the 
nature of the compact object. Indeed, observations in accretion 
disks orbiting around white dwarfs, neutron stars or black holes 
have shown a strong correlation between thek low and high 
frequencies QPOs (Mauche 120021 Psaltis, Belloni & van der 
Klis ll99^ . The relation is found to be same irrespective of the 
nature of the compact object. This strongly supports the idea 
of one and same general mechanism at hand in this accretion 
disks. The twin peaks ratio around 3:2 for the kHz-QPOs is nat- 
urally explained not only for black hole candidates or neutron 



stars but also for white dwarfs as reported by Warner, Woudt & 
Pretorius ( 2003 1. Indeed, the presence or the absence of a solid 
surface, a magnetic field or an event horizon play no relevant 
role in the production of the X-ray variability (Wiinands l2001> . 
The twin peak separation being either Av - v, for slow rotator 
(v. 300 Hz) or Av = v./2 for fast rotator (v. x 600 Hz) is 
also explained. Indeed, it was previously argued by Kluzniak, 
Lasota, Abramowicz & Warner (2005 1 that the 3:2 QPOs re- 
cently observed in white dwarf sources by Warner & Woudt 
(.2005,) have the same nature as the strong gravity 3:2 QPOs ob- 
served in neutron star and black hole sources. All of them are 
resonant accretion disk oscillations. Differences are attributed 
to different modes that are involved, to different mechanisms of 
resonance excitation, and to different modulations of the X-ray 
flux. 

To conclude, to date we know about 20 LMXBs containing 
a neutron star and all of them show kHz-QPOs. These QPOs 
can be explained by a mechanism similar to those exposed 
here. We need only to replace the gravitational perturbation by 
a magnetic one as described in Petri J2005> . However, in an ac- 
creting system in which the neutron star is an oblique rotator, 
we expect a perturbation in the magnetic field to the same order 
of magnitude than the unperturbed one. Therefore, the linear 
analysis developed in this paper has to be extended to oscilla- 
tions having non negligible amplitude compared to the station- 
ary state. Nonlinear oscillations therefore arise naturally in the 
magnetized accretion disk, leading to a shift in the resonance 
criteria and accounting for a change in the peak separation in 
relation with the accretion rate. 
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